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We perform extensive Monte Carlo simulations of a two-dimensional bent hard-needle model in both its chiral 
zig-zag and its achiral bow-shape configurations and present their phase diagrams. We find evidence for a 
variety of stable phases: isotropic, quasi-nematic, smectic-C, anti-ferromorphic smectic-A, and modulated- 
nematic. This last phase consists of layers formed by supramolecular arches. They create a periodic mod¬ 
ulation of the molecular polarity whose period is sensitively controlled by molecular geometry. We identify 
transition densities using correlation functions together with appropriately defined order parameters and 
compare them with predictions from Onsager theory. The contribution of the molecular excluded area to 
deviations from Onsager theory and simple liquid crystal phase morphology is discussed. We demonstrate the 
isotropic-quasi-nematic transition to be consistent with a Kosterlitz-Thouless disclination unbinding scenario. 

PACS numbers: 61.30.Cz, 64.70.M- 


I. INTRODUCTION 

In recent years, interest in two-dimensional and quasi- 
two-dimensional self-assembled structures in thin films 
has grown tremendously^^— . The possibility to fine- 
tune molecular ordering indeed makes thin films suit¬ 
able for myriad technological applications^, ranging from 
electronics^ and optics^ to biology^. The formation of 
liquid-crystal phases on two-dimensional surfaces is also 
key for various nanotechnological applications^^— . 

Control over self-assembly has already enabled 
the formation of two-dimensional aggregates with 
quasicrystalii^— , hexagonali^, crystali^, and liquid 
crystalii^— orders. Improving control on the sponta¬ 
neous formation of ordered structures, however, requires 
a deep understanding of how molecular geometry in¬ 
fluences supramolecular ordering. A common starting 
point for studying this relation are hard-core models, 
which recapitulate the static and dynamical properties 
of a wide range of phenomena, from the hard-sphere-like 
freezing of atomic liquids^S to quasicrystal formation^i. 
These models also reproduce the rich two-dimensional 
liquid crystal ordering of objects with high aspect ratios 
such as rods22i“— , rectangle a^^’^^ , spherocylinders2L2^, 
and ellipsoids22i2£. The good agreement between hard¬ 
core models and some experiments suggests that en¬ 
tropy alone can suffice to drive the formation of ordered 
structures^. 

High-aspect ratio molecules with a bent core22i2^ or 
a banana shap o^^’^^ assemble in an even richer set of 
morphologies, including biaxial nemati o^^i^^ , fan-shaped 
texture^^, and nematic with splay-bend deformation or 
conical twist-bend helix^. From the technological view¬ 
point these structures are quite interesting. Biaxial ne- 
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matic phases, for instance, exhibit ferro- or antiferro- 
electric properties'^. A plethora of models have thus 
been devised to understand the three-dimensional bulk 
behavior of these systems^^— . Yet relatively little theo¬ 
retical attention has been paid to related models in two 
dimensions'^— . 

Here, we consider the two-dimensional phase behavior 
of a bent-needle model in both its chiral zig-zag and achi¬ 
ral bow-shaped configurations (Fig. [T|). These two ver¬ 
sions display significantly different mesophases, for which 
we map out complete phase diagrams. Zig-zag molecules 
are known to assemble in either a nematic or a smectic-C 
phase depending on the packing density p^, and bow¬ 
shaped molecules have been found to display tetradic 
and nematic order—. Yet these studies did not clarify 
the role of topological defects and of thermal fluctua¬ 
tions on the long-range stability of the mesophases, which 
are of fundamental physical interest^^— . The relatively 
coarse sampling of configuration space previously used 
further left open the possibility that qualitative features 
of the molecular ordering may have been missed. We in¬ 
deed find that in addition to forming quasi-nematic and 
smectic phases, bow-shaped molecules present a stable 
modulated-nematic phase, as was previously predicted^^ 
and observed^2iS in three-dimensional systems. 

Our advances are not only made possible by the use 
of specialized Monte Carlo simulations on large systems, 
but also by the definition of appropriate order parame¬ 
ters and correlation functions. Our analysis thence ex¬ 
tends investigations of two-dimensional systems of hard 
needles^i^i^, hard spherocylinder o^^i^^ , and sphero- 
cylinders with a polar head^. These improved numerical 
results on two-dimensional liquid crystal formation are 
also compared to predictions from Onsager theory^. 

The paper is organized as follows. In Sec. [TTl we de¬ 
scribe the model and provide details on the Monte Carlo 
simulation procedure. Section m reviews the theory of 
quasi-long-range orientational order in two dimensions. 
In Sec. m we present the results of our analysis of the 
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FIG. 1. Bent-needle model molecule in chiral zig-zag (left) 
and achiral bow-shaped (right) configurations. 


simulation data. In Sec. El we summarize the Onsager 
Theory for two-dimensional hard-core objects and apply 
it to the bent-hard needle model. We conclude in Sec. 
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II. MODEL AND SIMULATION METHODS 
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A. Bent hard needles 

Our bent hard-needle model consists of a central line 
segment of length b to which two terminal line segments 
of equal length a are attached at a fixed angle a (Fig.[T]). 
The molecules can adopt a chiral zig-zag or an achiral 
bow-shape configuration, depending on whether a is de¬ 
fined on the same or on opposite sides of the central seg¬ 
ment. Note, however, that we only consider enantiomeri- 
cally pure systems in order to avoid chiral segregation^. 
Pairs of molecules interact via a hard-core exclusion po¬ 
tential but are infinitely thin, i.e., they can be infinitely 
close to one another but cannot overlap. Since the hard¬ 
core interaction potential is athermal, we set to unity 
the product of Boltzmann’s constant and temperature, 
= 1//3 = 1, without loss of generality. The total 
length of the molecule, L = 2a + b = 1, is used as the 
unit of length, also without loss of generality. Each model 
is thus completely determined by two parameters: a and 
a. In the following we let a vary from 0 to "k 12 (at a = 0 
both models are equivalent), but we fix a = 0.25 for zig¬ 
zag molecules and a = 0.35 for bow-shaped molecules. 
This choice maximizes the excluded area (to be defined 
below), and thus pushes down the isotropic-mesophase 
transition densities, making them computationally less 
costly to study. 

Figure [5] depicts the quasi-nematic and smectic phases 
for zig-zag and bow-shaped molecules as well as the 
modulated-nematic phase of bow-shaped molecules. This 
last phase, which consists of layers made up of arches 
formed by several molecules^^, is discussed in more de¬ 
tails in Sec. II V Cl 


FIG. 2. (colors online) Snapshots of the main mesophases 
identified in Monte Carlo simulations: (a) zig-zag molecules 
with a = 0.25, a = tt/S in both the quasi-nematic (p = 
14) and the smectic (p = 30) phases, and (b) bow-shaped 
molecules with a = 0.35, a = tt/IO, p = 20 in the quasi¬ 
nematic, with a — 0.35, a = tt/3, p = 30 in the smectic, and 
with a = 0.35, a = 7r/6, p = 26 in the modulated-nematic 
phases. The nematic director, the smectic surface normal 
and the modulated-nematic layers normal are indicated as 
red, blue, and green arrows, respectively. Note that in order 
to clearly illustrate the mesophase morphology, only a portion 
of the simulation box is displayed. 


B. Order Parameters 

We define the two-dimensional nematic order param¬ 
eter S by first introducing the tensor order parameter 

Q 

N 

Qo/3 = (iV-i^(2u>^-M), (1) 

i=l 

where is the cr-th Cartesian coordinate of the unit vec¬ 
tor pointing along the central segment of the i-th of N 
molecules and (...) denotes the ensemble average. The 
positive eigenvalue and corresponding eigenvector of Q 
give S and the nematic director n, respectively. The 
nematic order parameter S is canonically used to quan¬ 
tify the degree of molecular alignment along n - perfect 
alignment has 5" = 1. In two dimensions, however, only 
quasi-long-range orientational, i.e., quasi-nematic, order 
can exist (see Sec. IIII|) . and hence S' = 0 in the ther- 
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modynamic limit, N —>■ oo^. In order to characterize 
orientational order in the quasi-nematic phase, we resort 
to the orientational correlation function 

52 W = (cos[2(6»(0) - 6»(r))]), (2) 

where 6 is the angle between the central molecular seg¬ 
ment along u and a fixed axis, and r is the distance be¬ 
tween the centers of two molecules. The function 52 (^) 
thus monitors the spatial decay of orientational correla¬ 
tions. Note that because we work under periodic bound¬ 
ary conditions, correlations are radially truncated at half 
the edge length of the simulation box. 

A smectic liquid crystal can be thought of as a stack 
of parallel molecular layers of thickness d. It is there¬ 
fore possible to identify a density wave along the nor¬ 
mal to the layers. The smectic order parameter Asm is 
then the amplitude of this density wave. To determine 
Asm) we calculate the Fourier transform of the normalized 
density^. 


A(fe,d) 



N 

i=i 


( 3 ) 


where rj is the position of the j-th particle, and k = 
{2TT/d)k is the wave vector with wavelength d and unit 
vector k. To determine the actual distance d between 
smectic layers and the layer orientation fe, we calculate 
A{k,d) for different values of d and directions fc, and 
take its maximum Asm as the smectic order parameter—. 
Figure |3ja) illustrates this optimization procedure. 

Bow-shaped molecules are able to display polar order. 
To quantify it, we introduce the global polar order pa¬ 
rameter 



i=l 


( 4 ) 


where u* is a unit vector perpendicular to the molecular 
axis M* (see Fig.[T]). Since Pg cannot distinguish between 
isotropic and anti-polar order - in both cases Pg = 0 
-, other order parameters are needed to characterize an 
anti-polar state. In Sec. IIV Cl we describe the modulated- 
nematic phased that bow-shaped molecules form at in¬ 
termediate a. In this mesophase molecules form arches in 
which the orientation of the molecules’ central segments 
- and thus also the polar vector v (see Fig. [1]) - display 
a periodic modulation along one particular spatial direc¬ 
tion. To quantifiy this spatial modulation, we introduce 
its amplitude as an order parameter Tm, which we de¬ 
termine in full analogy to the definition of the smectic 
order parameter Asm- We define the spatially dependent 
polarity, Pi(t’) = ~ calculate the 

amplitude of its Fourier transform 
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FIG. 3. (colors online) Illustration of the maximization proce- 
dnre for determining the smectic (a) and modulated-nematic 
(b) order parameter, (a) Fourier transform of the density, 
A(fe,d), plotted versus orientation angle (pi and wavelength 
d of the wave vector k = (27r/d)fc with fc = [cos<?ii,sin0i]. 
The data refer to zig-zag molecules with a = 0.25, a = 7 r/ 2 , 
and p = 30. The maximization gives Asm = 0.83. (b) Fourier 
transform of the polar order parameter, ^(7i, A), plotted ver¬ 
sus orientation angle (j )2 and wavelength A of the wave vector 
h = (27r/A)/r with h = [cos<?!> 2 , sinThe data refer to bow¬ 
shaped molecules with a = 0.35, a = 7r/13, p = 24. The 
maximization gives ^'m = 0.62. In both cases the procedure 
gives rise to a clear maximum, which is used for further anal¬ 
ysis. Note that only a portion of the domain explored during 
the optimization procedure is displayed here. 


where h = {2 t:/\) h is the wave vector with wavelength 
A and unit vector h. Note since 9 + tt/2 quantifies the 
orientation oft),'!' describes periodic variations in v. We 
evaluate 'I'(/i, A) for different values of A and directions h 
and take the maximum Ifm as the order parameter for the 
modulated-nematic phase. The optimization procedure 
is illustrated in Fig. [3Kb). 

In order to have additional information on the struc¬ 
tural properties of the modulated-nematic phase, we 
monitor polar correlations along the optimal direction 
hni using the polar correlation function 


N 

7=1 


Ah', 


( 5 ) 
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Here, = |rj — rj \ ■ is the projection of the vector 
joining the centers of the i-th and j-th molecules onto 
the direction of /im- 

C. Monte Carlo simulations 

In order to study the phase behavior of the zig-zag and 
bow-shaped models, we perform Monte Carlo simulations 
mostly in the constant-but also in the constant- 
NPT ensemble under periodic boundary conditions^. 
Most simulations have N = 2000 in a square box of area 

Constant-simulations are performed at hxed re¬ 
duced density p = N/V and consist of up to 3.0 x 10^ 
sweeps, where a Monte Carlo sweep comprises N in¬ 
dependent trial displacements. The basic Monte Carlo 
move in the NVT ensemble consists of either translat¬ 
ing or rotating (randomly chosen with equal probabil¬ 
ity) a randomly chosen molecule. The maximum at¬ 
tempted displacement is preliminarily adjusted in such 
a way that an acceptance rate of approximately 50% is 
achieved, but during the production runs the maximum 
step sizes are hxed - dynamical adjustments would vi¬ 
olate detailed balance^. The hard-core nature of the 
pair interaction simplihes the Metropolis acceptance rule: 
the displacement of a molecule is accepted if it does not 
generate an overlap, and is rejected otherwise. In or¬ 
der to check that the sampled conhgurations correspond 
to equilibrium (and not metastable) states, different ini¬ 
tial conhgurations - including isotropic, perfect polar and 
anti-polar nematic, as well as perfect smectic states - 
are used. Note that reaching equilibrium and properly 
sampling states sometimes require up to several million 
Monte Carlo sweeps, leading to computational times of 
up to four weeks on a single-thread of an Intel Xeon 
X5550 machine with a 2.66 GHz CPU. The minimum 
simulation time was about one week. 

Constant-A^PT simulations also include changes in V 
in order to keep the system pressure P constants^. A 
Monte Carlo sweep then comprises an average of N inde¬ 
pendent single-molecule trial displacements and one trial 
volume change. In a trial volume change, V is modihed 
by an amount A In U and all the molecule positions are 
rescaled accordingly. If no overlap occurs in the resulting 
conhguration, the move is accepted with probability 

p{V V + AV) = exp[-PAU-h(Ar-^I)ln(I-HAU/U)], 

(7) 

and is rejected otherwise. 

For the bow-shaped molecules with intermediate a, 
constant-AUT simulations either result in ferromorphic 
or in anti-ferromorphic states, depending on the initial 
conhguration. In order to identify the equilibrium ground 
state, we thus also performed a slow pressure anneal¬ 
ing. This process was achieved through a sequence of 
constant-APT simulations starting from the isotropic 
regime and increasing P in steps of 10-20%. At each step 


p was equilibrated for at least 2.0 x 10® sweeps. The op¬ 
posite procedure was used to check for hysteresis. This 
study revealed that the anti-ferromorphic phase is the 
equilibrium state, as we discuss in Sec. II V Cl 

Since the interaction between bent needles is short- 
ranged, we use a combination of linked and Verlet lists 
to check for the overlap of molecules after a Monte Carlo 
trial. This speeds up the simulation s®®:®^ . The Verlet list 
is built by wrapping a spherocylinder body around each 
molecule. This non-trivial shape is obtained by gluing 
together three spherocylinders, one for each molecular 
segment. The Verlet list is then filled with the neighbor¬ 
ing molecules whose bent spherocylinders overlap with 
the one under consideration. The bent spherocylinder of 
a molecule remains fixed for many Monte Carlo steps, 
but has to be reconstructed when a rotational or a trans¬ 
lational Monte Carlo move brings the molecule out of its 
bent spherocylinder. Shrinking the spherocylinder radius 
for the Verlet list considerably reduces the time needed to 
detect potential overlaps, but increases the rate at which 
the lists need to be updated^. Before starting a Monte 
Carlo simulation, we thus first determine the spherocylin- 
der radius that optimizes algorithmic performance. 

The linked list is built by using a square decomposition 
of the simulation box. The simulation box is divided 
into smaller sub-boxes of a side length roughly given by 
the particle length plus the spherocylinder radius of the 
Verlet list. The Verlet list in a given sub-box can then be 
built by only considering molecules within that same sub¬ 
box as well as within the eight neighboring sub-boxes. 

D. Cluster moves 

In the vicinity of the isotropic-quasi-nematic tran¬ 
sition, small clusters of very close and well aligned 
molecules develop in the isotropic phase. Once these rela¬ 
tively small but highly packed clusters of particles form, 
molecule orientations can get kinetically locked. Equi¬ 
libration, however, requires overcoming the high free- 
energy barriers associated with aligning these clusters 
which considerably slows down equilibration. A way to 
alleviate this problem relies on the fact that Monte Carlo 
simulations need not be tied to a local and thus physical 
dynamics. Instead of moving particles one at a time, one 
can define collective moves that identify groups of corre¬ 
lated particles and then move them as a single objec t®®’®"^ . 
However, in order to satisfy detailed balance and ensure 
that reverse displacements are also possible, the opera¬ 
tion must be done in a probabilistic way. 

Here, we build clusters of particles by introducing an 
artificial attractive potential Uf^Cij) that links particles 
togethe r®®!®® . The potential form is formally arbitrary, 
but is most successful if it captures the nature of corre¬ 
lations within a cluster. Our choice is 

w/(eb) = -eH(A6*i„ax - A6'y )H(Armax - Ar^), (8) 

where H is the Heaviside step function, and Armax and 
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FIG. 4. Illustration of the calculation of the excluded area 
for bow-shaped molecules with a = 57.3°, a = 0.35 and the 
same orientation. The excluded area is obtained by repeat¬ 
edly inserting the second molecule inside the box and checking 
for overlaps. The shaded region is obtained by placing a dot 
at the position of the center of mass of a randomly inserted 
particle if it overlaps with the fixed particle. 


A^max are tunable thresholds. Typical values are e = 1.0, 
Ar^ax = 0.3 and A6*i„ax = 6.0°. For a given configura¬ 
tion, the linking probability between particles i and j 
is 

Pij = niax[0,1 - exp(/3u/(eij))]. (9) 

Molecules that are close and well aligned are thus linked 
with high probability, whereas molecules outside the 
range of Armax and A61max are not linked at all. A clus¬ 
ter is built by choosing a particle i at random and by 
attempting to build links with its neighbors. If a link 
between i and j is formed, j becomes a member of the 
cluster formed around molecule i. Once the cluster has 
been formed, we perform a trial move by translating or 
by rotating it around its center of mass. If the trial move 
does not generate any overlap with other molecules, it is 
accepted with probability 

W (old new) = min[l, exp(/3(t7new - C^oid))], (10) 

where Ui, = '^i is the artificial interaction 

energy between the cluster and its environment in either 
the V =old or new state. The interface is defined by 
all the particles i inside and the particles j outside the 
cluster that contribute to Uy. 

E. Evaluation of excluded areas 

In Sec. 0we use Onsager theor y^^'^^ as an alternative 
method for calculating phase diagrams. This approach 
heavily depends on the the concept of excluded area. 
Given two molecules with a fixed relative orientation A0, 
then the excluded area Aexi(A6>) is defined as the portion 
of the plane surrounding a molecule that is not accessible 


to the other. If the center of the second molecule is inside 
the area excluded by the first, then molecules overlap. 
Figure m shows an example of the shape of the excluded 
area for two bow-shaped molecules with the same orien¬ 
tation (A6> = 0). For hard-core potentials the excluded 
area is proportional to the second virial coefficient in an 
expansion of the free energy in powers of density. Hence, 
Aexi provides a microscopic description of the balance 
between positional and orientational entropy. 

We evaluate the excluded area by Monte Carlo inte¬ 
gration. One molecule is fixed at the center of a box of 
area Abox- The size of the box is chosen such that it 
is impossible for a molecule outside the box to overlap 
with the fixed molecule. Another molecule with fixed 
relative orientation A0 is then inserted at random in the 
box. The process is repeated Atrial times, keeping track 
of the number of overlaps Aoveriap that occurred dur¬ 
ing the whole process. In the end, the excluded area is 

Aexl(A^) = Abox (Aoverlap/Atrial)- 


III. QUASI-LONG-RANGE ORIENTATIONAL ORDER 

According to the Mermin-Wagner theorem the spon¬ 
taneous symmetry breaking of a continuous order pa¬ 
rameter is always suppressed by fluctuations in dimen¬ 
sions d < 2 for systems with sufficiently short-ranged 
interactions^. Nevertheless, as first pointed out by 
Kosterlitz and Thouless (KT), a phase transition in d < 2 
is still possible between, on the one hand, a disordered 
phase wherein the correlation function of the order pa¬ 
rameter decays exponentially and thereby only exhibits 
short-range order, and, on the other hand, a phase with 
quasi-long-range order wherein the correlation function 
decays as a power law in distance r— . 

Since it will be important for our analysis in Sec. II V A1 
we summarize the main KT results. The description 
starts from the free energy associated with distortions 
of a molecular orientational field, 

F = ^ J [ye{r)r d\, ( 11 ) 

where the angle 9{r) measures orientation of the molecule 
with respect to a fixed axis, and K is the Frank elas¬ 
tic constant. In a more general theory, Eq. (ED should 
include two different elastic constants, one for splay 
and one for bend deformations, but on sufficiently large 
length scales they renormalize to the same value^i. The 
orientational correlation function ( 72 (^) introduced in Eq. 
([2]) quantifies the observed orientational order. Accord¬ 
ing to Kosterlitz and Thouless, quasi-long-range orienta¬ 
tional order (here, a quasi-nematic phase), results from 
the competition between the free energy needed to cre¬ 
ate topological defects and the entropy gained when these 
defects unbind and are thus free to move^^. Disclination 
unbinding takes place at the critical value of the Frank 
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elastic constant^^, 


-kKc 

SksT 


( 12 ) 


We stress that Kc is a scale-free quantity that locates the 
transition between short-range and quasi-long-range ne¬ 
matic order in the thermodynamic limit^^. The critical 
value is the result of a balance between the disclination 
energy and entropy, and in two dimensions both have the 
same logarithmic dependence on system size. If if < Kc, 
isolated disclinations are found, which leads to an overall 
isotropic state characterized by an exponential decay of 
the orientational correlation function. If if > Kc, discli¬ 
nations are bound in pairs and the orientational corre¬ 
lation function is thus expected to decay algebraically, 
52 (^') oc r~^, with an exponent^S 


7i = 2kBT/TrK. (13) 


In the limit a = 0, (or for a = 0) both configurations of 
the bent-needle model reduce to a straight needle, which 
is known to undergo an isotropic-quasi-nematic transi¬ 
tion via disclination unbinding^i^. We demonstrate in 
Sec. IIV Al that the isotropic-quasi-nematic transition in 
the bent-needle model is also consistent with the KT pic¬ 
ture. 



IV. RESULTS 

In this section we present the results of the Monte 
Carlo simulations described in Sec. im We detail how 
we identify the isotropic-quasi-nematic transition as well 
as the subsequent quasi-nematic-smectic transition. We 
also introduce the modulated-nematic phase and summa¬ 
rize our results in phase diagrams. 

A. Isotropic-quasi-nematic transition 


FIG. 5. (colors online) (a) Radial dependence of the orienta¬ 
tional correlation function for several densities for (a) zig-zag 
molecules with a = tt/S and a = 0.25, and (b) bow-shaped 
molecules with a = 27r/5 and a = 0.35. Dashed lines in (b) 
show the fits to an exponential function for p = 15 and 16 and 
to 52 (r) oc r~^ for p = 18, 20, and 23. Fit exponents are given 
in Fig. [6] The inset in (a) shows a portion of the simulation 
box at p = 11 , close to the isotropic-quasi-nematic transition. 
The inset in (b) details the local packing of molecules, which 
results in oscillations of 52 (r). 


Upon increasing density, both zig-zag and bow-shaped 
molecules form a quasi-nematic phase with quasi-long- 
range orientational order . The quasi-nematic phase and 
the isotropic-quasi-nematic phase transition via disclina¬ 
tion unbinding are evidenced by the transition from an 
exponential decay of 52 (j’) at low p to a power-law de¬ 
cay, g2{'i") oc r~^, as p increases (see Fig. El). As per 
the discussion in Sec. mi the quasi-nematic phase is ex¬ 
pected to be stable against spontaneous disclination un¬ 
binding when -kK/ {SksT) > 1. We obtain the Frank 
elastic constants from the relation rj = 2A:bT juK, where 
the power-law decay of 52 (*") is obtained from a linear 
fit of \n.g2{r) versus Inr over the range 0.5 < r < 4.0. 
Note that the lower threshold is necessary because 52 {r) 
deviates from the power law at small r, while the higher 
threshold is chosen so as to exclude correlations resulting 
from the use of periodic boundary conditions. Note also 
that bow-shaped molecules display correlation functions 


with fairly large oscillations as a result of the local pack¬ 
ing structures, which do not exist for straight needles [see 
inset in Fig. EKb)]. 

The results for the reduced Frank elastic constant 
ttK/ i^ksT) are plotted in Fig.Elfor both molecule types. 
Since the results scale nearly linearly with p, we iden¬ 
tify the isotropic-quasi-nematic transition densities piN 
as the intersection between linear fits to data the points 
and irKc/iSksT) = I—. 

For zig-zag molecules, the quasi-nematic phase is sys¬ 
tematically destabilized by bending the terminal seg¬ 
ments [Fig. ini[a)]. As the central- to-tail angle a in¬ 
creases, the transition is thus pushed to higher densities. 
Bow-shaped molecules, however, show a non-monotonic 
trend of the transition density piN with a [Fig. EKb)]. 
For small a the transition density increases with a, 
while for a > tt!21 the trend is inverted. In the range 
7r/I3 ^ a < 7r/3 the quasi-nematic phase is unstable 
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(a) 




FIG. 6. (colors online) Density evolution of the reduced Frank 
elastic constant K obtained from the power-law decay of g 2 (r) 
for (a) zig-zag molecules with a = 0.25 and (b) bow-shaped 
molecules with a = 0.35. The isotropic-quasi-nematic transi¬ 
tion densities, piN are obtained from the intersections of the 
linear fits to the numerical results with the line nK/SksT = 1. 
Where not shown, errorbars are smaller than the symbols. 


with respect to the modulated-nematic phase, and the 
power-law scaling analysis is then inapplicable (see Secs. 
HVCl and UVDll . 

We try to understand the difference in the two molec¬ 
ular geometries by examining the size Aexi of the ex¬ 
cluded area for perfectly parallel molecules (see Fig. El). 
Like Pin, it monotonically increases with a for zig-zag 
molecules, but reaches a maximum at a w 7r/3 and then 
decreases for bow-shaped molecules. Large values of Aexi 
suggest that parallel molecules have to pack locally to 
come close to one another, as evidenced by the undu¬ 
lations of g 2 {'i’) in Fig. [5Kb). This packing constraints 
translational freedom, and thus reduces the translational 
contribution to the entropy. Because the translational 
contribution normally compensates for the loss of ori¬ 
entational freedom in the nematic phase, higher p than 
usual are needed for this effect to be significant, and as 
a result piN increases. 



FIG. 7. Excluded area versus central-to-tail angle a for par¬ 
allel zig-zag (continuous line) and bow-shaped (dashed line) 
molecules. 


B. Quasi-nematic-smectic transition 

Smectic order consists of a periodic arrangement of 
particle positions along one direction, which leads to a 
well-defined density wave along the corresponding wave 
vector. Within smectic layers, however, translational 
order is absent. Although long-range translational or¬ 
der is not expected in the thermodynamic limit of two- 
dimensional systems, smectic order is obervable on suf¬ 
ficiently small length scales. We get back to this point 
below. For now, we consider the smectic order parameter 
defined in Eq. dS]) as a function of density for both the 
zig-zag and bow-shaped molecules in Fig. [51 The Monte 
Carlo simulation data are fitted with 

/(p) = 1/2 -b arctan[/i(p - pnsKI/tt, (14) 

where pNS and h are fit parameters. Hence, the quasi¬ 
nematic-smectic transition density pNS is defined as the 
point of maximum slope of /(p). The smectic order pa¬ 
rameter for bow-shaped molecules does not assume large 
values even at the highest densities simulated, which is 
likely the result of out-of-layer fluctuations and of the in¬ 
stability of two-dimensional smectic order described. Fits 
to the simulation results nonetheless provide an estimate 
for PNS- Note that the values of the transition densities 
identified in this way are close to the highest values of 
the densities explored in our simulation. 

As could be seen in Fig. El and was previously noted 
in Ref. IH, zig-zag molecules form a smectic-C phase, 
wherein the nematic director is tilted with respect to 
the layer normal. Bow-shaped molecules instead arrange 
in an anti-ferromorphic smectic-A structure, wherein the 
polar vector v adopts an opposite orientation in adja¬ 
cent layers. This unusual ordering can be rationalized by 
extending a packing argument developed for the anti- 
ferrolectric ordering of V-shaped molecules^, i.e., for 
b = OiS. Because the excluded area of two molecules 
is smaller in the anti-parallel (D^*^ = than in 
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FIG. 8. (colors online) Smectic order parameter Asm plotted 
versus density p at several central-to-tail angles a for (a) zig¬ 
zag molecules with a = 0.25 and (b) bow-shaped molecules 
with a — 0.35. Data points are obtained from Monte Carlo 
simulations. Lines are a fit of Eq. d to the data using the 
transition density pNS and h as fitting parameters (see phase 
diagrams of Fig. [T^] for results on pns). 


the parallel arrangement, the former is en- 

tropically favored. Molecules can also more easily pen¬ 
etrate into neighboring layers in the anti-ferromorphic 
than in the ferromorphic smectic phase. This effect en¬ 
hances out-of-layer fluctuations and thus entropically fa¬ 
vors the anti-ferromorphic smectic ordering as well. This 
behavior is similar to that of three-dimensional bent-core 
molecules^^. 

No significant smectic ordering is found for values of a 
smaller than those given in Figs.lSKa) and (b). For zig-zag 
molecules at even smaller a, we extrapolate the quasi¬ 
nematic-smectic transition to take place at pNS that are 
inaccessible within a reasonable computational time; for 
bow-shaped molecules, however, the existence of a smec¬ 
tic phase at q; < tt/S, even for large p, is unclear. 

As mentioned above and discussed in detail by Toner 
and Nelson in Ref. [tI, two-dimensional long-range smec¬ 
tic order should not be thermodynamically stable. One 
expects instead the following scenario. At a given p trans¬ 
lational order is only disturbed by phonon fluctuations 


in regions with linear dimension smaller than a charac¬ 
teristic length which is the mean distance between 
thermally induced dislocations. These dislocations thus 
destroy translational order on length scales larger than 
but correlations in the layer orientations persist and 
exhibit a long-range algebraic decay. Although this sce¬ 
nario is physically reasonable and may explain the weak 
smectic ordering of bow-shaped molecules, we were un¬ 
able to test it against our simulation data because of the 
limited range of computationally accessible N. Instead, 
we And smectic order to be fully stabilized in our simu¬ 
lations. 


C. Modulated-nematic phase 

Although at small a and for a > tt/S bow-shaped 
molecules form a quasi-nematic phase [Fig. [2][b)], for 
intermediate a, they also equilibrate in a modulated- 
nematic phase. In this phase no overall orientational 
or positional order exists [Fig. |9l[a)]. Instead, it shows a 
different kind of supramolecular arrangement, wherein 
the orientation of a series of molecules varies gradu¬ 
ally along arches that form approximate half-circles, and 
these arches themselves form layers. Periodic order ex¬ 
ists along the layer normal but the mean polar vector (v) 
in one layer is antiparallel to that of a neighboring layer. 
This arrangement destabilizes the quasi-nematic phase 
and its algebraic orientational order. As shown in Fig. 
ini[b), the occurrence of the modulated-nematic phase can 
even result in a reentrant exponential decay of g 2 (r) with 
increasing p. For 7r/13 ^ a tt/S, however, we found 
no evidence for a power-law decay of g 2 (f) in the whole 
range of explored densities and a direct transition from 
the isotropic to modulated-nematic phase occurs. 

In Fig. llOf a) the modulated-nematic order parameter 
defined in Eq. ([5]) is shown as a function of density 
for several values of a. In analogy with the treatment 
for the identification of smectic order in Sec. lIVBi we fit 
Iflni (as obtained from the Monte Carlo simulation data) 
with a trial function similar to Eq. da, 

f{p) = 1/2 -b arctan[;(p - p^)]/tt (15) 

where I and pm are fit parameters. In order to clarify the 
structural properties of the modulated-nematic phase, we 
consider the polar correlation function ^^(r) as defined 
in Eq. Because of the periodic modulation in the 
molecular polar vector v along the layer normal, Pi (r) 
becomes a periodic function when the layered structure 
is well established. In particular, since molecules in ad¬ 
jacent layers have opposite polarization, p['(r) shows a 
minimum at a distance corresponding to the layer thick¬ 
ness followed by a maximum, which results from correla¬ 
tions with the next-nearest-neighbor layer. In Fig. ITUI b) 
we show Pi(r) for a = tt/IO and several values of p. 
As expected, the periodic modulation of g^l (r) becomes 
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FIG. 9. (a) Snapshot of bow-shaped molecules in the 

modulated-nematic phase for a = n /8, a — 0.35, and p = 20. 
The complete simulation box with N = 2000 molecules is 
shown. The red arrows indicate the polar vector v perpendic¬ 
ular to the molecular orientation in the center of the layers. 
The blue arrow indicates the direction of the optimal wave 
vector hm- The radius of the blue circle is half the opti¬ 
mal wavelength Am as defined in Eq. (0 and obtained from 
the maximization procedure illustrated in Fig. [Sj^b). For this 
configuration « 0.25. (b) The orientational correlation 
function <72 (r) for bow-shaped molecules with a = tt/IS de¬ 
cays exponentially at p = 6 , algebraically at p = 11 , and 
exponentially again at p = 13 due to the appearance of the 
modulated-nematic phase. 


stronger with increasing density indicating the progres¬ 
sive development of layers. This behavior is also observed 
for all other values of a. 

Our Monte Carlo data show that modulations in (r) 
are already present at relatively small p, whereas the or¬ 
der parameter has not yet increased much. For in¬ 
stance, for a = 7r/10 it can be seen in Fig. [TUT bl that 
gi{r) is already weakly modulated at p = 14, while the 
corresponding order parameter is only ~ 0.2. Simi- 



FIG. 10. (colors online) (a) Modulated-nematic order param¬ 
eter i&in plotted versus p for bow-shaped molecules at several 
values of a. The Monte Carlo simulation results are fitted 
with Eq. (1151) . From the fit we obtain pm = 29.95, 19.97, 
18.71, 21.42, 23.38, and 45.57 for a = 7r/25, 7r/16, 7r/13, rr/lO, 
7 r/ 8 , rr/b, and tt/4, respectively, (b) Polar correlation func¬ 
tion Qiir) for a = tt/IO and several values of p. As density 
increases, the amplitude of the modulation of increases, 
which we take to be a signature of the modulated-nematic 
phase. 


larly, we find 'km Ri 0.25 for the configuration shown in 
Fig. UKa), even though a remarkable degree of layering is 
already clearly established. In other words, the layered 
structure of the modulated-nematic phase is established 
continuously, which makes it difficult to clearly establish 
a transition density. For this reason, we approximate the 
transition density pm into the modulated-nematic phase 
by using the threshold fkin = 0.2 on the order parame¬ 
ter. The same threshold is used for all a. This particular 
value is chosen in such a way that pm matches reason¬ 
ably well the values of the density at which we observe 
the reentrant exponential decay of P 2 (?') at a < tt/IS. 

The maximization procedure illustrated in Fig. EKb) 
for determining the order parameter 'Flu also provides 
the typical distance between two layers in the modulated- 
nematic phase, i.e., Ani/2 defined in Eq. (0). In Fig. [TT] 
we show how this distance evolves with a and com¬ 
pare it with the projected length of the molecule along 
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a. 

FIG. 11. (colors online) Distance between two adjacent lay¬ 
ers, Am/2, plotted versns a (points). The results for Am are 
obtained at p > pm, bnt are observed to be fairly insensi¬ 
tive to changes in p [inset (3)]. The continnous red line is 
the projected length of the molecnle onto a direction along 
the central molecular segment. The vertical line at q = 7r/4 
approximately divides the plot in two regions. The first re¬ 
gion with a < 7r/4 has A/2 > Lc and molecules are ar¬ 
ranged in the supra-molecular structure of the modulated- 
nematic phase illustrated in inset (1). The second region with 
a > 7r/4 has A/2 < Lc, where molecules locally arrange with 
anti-ferromorphic order with adjacent layers slightly inter¬ 
penetrated, as illustrated in inset (2). The dotted line is a 
guide to the eyes. 


the direction parallel to the central molecular segment, 
Lc = b + 2acos(a). At small a, Am > Lc, which cap¬ 
tures the supra-molecular structure of the modulated ne¬ 
matic phase. As a increases, however. Am decreases, and 
it becomes comparable to Lc at a ~ 7r/4, indicating 
that molecules locally arrange in anti-ferromorphic or¬ 
der, as demonstrated by the inset (2) of Fig. [Tl] which 
is then favored over the supra-molecular ordering of the 
modulated-nematic phase. For a > Tr/4 the layering 
distance becomes smaller than Lc, indicating that ad¬ 
jacent anti-ferromorphic layers on average slightly inter¬ 
penetrate, as discussed in Sec. IIVBI 

The occurrence of a three-dimensional, spontaneously 
formed, modulated-nematic phase for banana-shaped 
mesogens and of the two-dimensional modulated-nematic 
phase for bow-shaped molecules, likely results from 
“pathological elasticity”—, which has been predicted to 
be a consequence of the molecular curvature radius. In 
the standard Frank elastic theory, splay, bend, and twist 
elastic constants are indeed assumed to be positive in 
order to ensure a ground state with uniform nematic 
order. A non-uniform nematic ground state can thus 
be explained as resulting from a negative elastic con¬ 
stant, e.g., the bend constant, which is not forbidden 
by symmetry. Fourth-order terms in the elastic free en¬ 
ergy are then needed to stabilize the modulated phased. 


Our results thus indicate that there is an upper limit for 
the molecular curvature, corresponding in our model to 
a ~ 7r/4, beyond which the supra-molecular structure of 
the modulated-nematic phase becomes unstable in favor 
of the anti-ferromorphic arrangement. 

We already discussed in Sect. IIV 51 that long-range or¬ 
der in two-dimensional smectics is not stable. In Fig. 
[ni[b) in the center a dislocation in the layering of the 
supramolecular arches is visible. This might be an indi¬ 
cation how, in analogy to two-dimensional smectics, also 
the layered structure is destabilized by the proliferation 
of dislocations in sufficiently large systems. 


D. Phase diagrams 

In Fig. [T^ we show the simulated phase diagrams for 
both the zig-zag and the bow-shaped molecules obtained 
using the approaches described in Secs. IIV A[ IIVBI and 
IIV Cl The dashed red and dotted blue lines show the pre¬ 
dictions from Onsager theory for the transition densities 
Pin and pNS, respectively. We will review Onsager theory 
in Sec. 0 

According to the phase diagram in Fig. fT^ a). termi¬ 
nal segments of zig-zag particles destabilize the quasi¬ 
nematic phase but stabilize the smectic phase. This 
qualitative trend is captured by Onsager theory. Our 
simulation results are in rough agreement with previ¬ 
ous studies^i^ but also exhibit some significant discrep¬ 
ancies. In particular, in our phase diagram both the 
isotropic-quasi-nematic and the quasi-nematic-smectic 
transitions appear at higher p (Fig. IT^ a)). In addition, 
for a ss 7r/2, the isotropic and smectic phases are well 
separated by the quasi-nematic phase, while the latter 
was not observed at all in previous studies. Most likely, 
the lower transition densities reported in Refs. andlSTi 
follow from using the nematic order parameter for detect¬ 
ing the isotropic-quasi-nematic transition and from con¬ 
sidering relatively small system sizes. As already men¬ 
tioned, the nematic order parameter in a quasi-nematic 
phase decreases with system size^. 

Bow-shaped molecules exhibit a remarkably rich be¬ 
havior. 

(1) For small a, the low-density behavior is similar to 
that of zig-zag molecules. The transition density piN (red 
circles) increases with increasing a and is underestimated 
by Onsager theory. Quasi-nematic order in this small-a 
region is also destabilized by further increasing the pack¬ 
ing density. Instead of the emergence of a smectic phase, 
however, a modulated-nematic phase takes precedence. 
We did not specifically explore the phase behavior of 
bow-shaped molecules with a < 7r/25, but we expect 
the modulated-nematic phase to appear at ever higher 
densities as a decreases, as indicated by the dotted green 
line. For symmetry reasons the bent-nematic phase can¬ 
not exist for a = 0. 

(2) For 7r/13 < a < 7r/4, both the quasi-nematic and the 
smectic orders are destabilized by the modulated-nematic 
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(a) 




a 

FIG. 12. (color online) Phase diagram for (a) zig-zag 
molecules with a = 0.25 and (b) bow-shaped molecules with 
a — 0.35. Points represent the isotropic-quasi-nematic tran¬ 
sition densities (circles), quasi-nematic-smectic transition 
densities (diamonds) and the quasi-nematic- or isotropic- 
modulated-nematic transition (squares) as identified from 
the procedures described in Secs. IIV A1 II V Bl and IIV Cl re¬ 
spectively. The gray area in (b) marks the transition be¬ 
tween modulated-nematic and quasi-nematic phases in a re¬ 
gion where none of the identified mesophases is found to be 
stabilized. Continuous lines are guides to the eyes. Where 
not shown, errorbars are smaller than the marker size. Dashed 
and dotted lines represent, respectively, the isotropic-nematic 
and the nematic-smectic transition lines predicted by Onsager 
theory (Sec.|V|. 


phase at all p explored. The direct isotropic-modulated- 
nematic transition appears at increasing densities with 
a. This phase is simply not captured by Onsager theory 
as formulated in Sect. El 

(3) For 7r/4 < a < tt/S, which is indicated in the phase 
diagram as a gray area, we find no clear evidence for any 
of the mesophases known to be formed by bow-shaped 
molecules. In this region the distance between two adja¬ 
cent layers approaches the projection of the total molecu¬ 
lar length along the molecular central segment, resulting 


in a strong competition between the formation of supra- 
molecular layers and anti-ferromorphic domains. The 
investigation of configurations at densities much higher 
than the ones we could reach would be necessary to clar¬ 
ify what mesophase, if any, is stable in that system. 

(4) For 7r/3 < a < 7r/2, the molecular curvature is 
too large to induce a spontaneous bending of the ne¬ 
matic director and the corresponding layer formation. 
Quasi-nematic order then becomes stable again, but this 
time the transition densities decrease with increasing a. 
The isotropic-quasi-nematic transition is highly overes¬ 
timated by Onsager theory. The theory strongly relies 
on the excluded area, as discussed in Sec. IIIEI but the 
excluded area of bow-shaped molecules is minimal in the 
anti-parallel configuration. In our simulations we instead 
observe clusters of molecules packed in parallel. Because 
this arrangement involves structural correlations between 
more than two molecules, it is not surprising that On¬ 
sager theory, which only takes into account two-particle 
correlations dramatically fails. 

Note that we did not study molecules with a = -k 12. 
With this central-to-tail angle parallel molecules cannot 
anymore be shifted into each other to form closely packed 
clusters. A different phase is thus expected. Tetradic 
order has indeed been observed in related models^. 


V. ONSAGER THEORY FOR PHASE DIAGRAMS 


In Ref. [13, a density functional theory in mean-field 
approximation was used to determine the phase behav¬ 
ior of zig-zag molecules. Based on the second-order virial 
expansion of the free energy introduced by Onsager, pre¬ 
dictions for both Pin and pNS were made^. In the follow¬ 
ing, we summarize and extend the results of this Onsager 
theory. 

We start with the free energy functional 


PF = J dr J dijj{\a.p{r,uj) — 1} 

J J diPip{fi,iPi) 

^ I I ‘^‘^2P(^2,W2)/M(fb2,Wi,W2) . (16) 


where p(r, w) is the local number density in terms of par¬ 
ticle position r and particle orientation u). The Mayer 
function /m of the pair potential is simple for hard-core 
particles. It is zero except when particles overlap, where 
/m = -1- 

In the isotropic and nematic phases, positional order 
does not exist and one has p(r, w) = pf{uj). In the second 
term of the Onsager free energy functional dUD, this in¬ 
troduces the excluded area Aexi(Q;) in terms of the angle 
a between the two molecular orientations in two dimen¬ 
sions. Expanding orientational distribution /(w) = f{9) 
and Aexi(cr) into Fourier modes and minimizing with re¬ 
spect to the Fourier amplitudes of f{9), gives a series of 
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possible bifurcation densities 


(n) ^ 

Pb 


(17) 


Here, A„ is the n-th Fourier amplitude of the excluded 
area Hexi(a), 


A 


n 


2m 

TT 


7r/m 

J Aexi(o) COS (mna) da . 
0 


(18) 


We have introduced an index m = 2 for zig-zag molecules 
and m = 1 for bow-shaped molecules because their re¬ 
spective excluded areas either have a period tt or 2tt. 
We determine the Fourier coefficients numerically 
by means of the Monte Carlo integration technique de¬ 
scribed in Sec. Ill El We choose the lowest positive value 
of as the transition density piN, which is realized at 
n = 1 for zig-zag molecules and at n = 2 for bow-shaped 
molecules. The resulting transition lines are plotted as 
red dashed lines in the phase diagrams of Fig. [TH 

To determine the nematic-smectic transition density 
PnS) we choose a periodic modulation for the density 
along the z axis, p{r) = p{z) = p{z + d), and assume per¬ 
fect alignment of the molecules as in Ref. [13, because the 
nematic order in the smectic phase is typically very high. 
The evaluation of the free energy functional now involves 
an excluded distance de-^\{z, 0), where 9 is the orientation 
angle of the central molecular segment with respect to the 
smectic layer normal. The excluded distance is related to 
the excluded area by = f dz dex\(z,9). Taking into 
account only the first Fourier mode of the density mod¬ 
ulation in the free energy and minimizing with respect 
to the Fourier amplitude, we obtain an equation for the 
nematic-smectic transition density pnSi 


1 + Pns 


COs{2TTz/d)de:^l{z, 9) dz = 0 , 


(19) 


where d is the smectic period. The tilt angle 9 and the 
period d are also determined by minimizing the free en¬ 
ergy functional. For detai ls of the cal culation, we refer 
the reader to the work of IVarga et al.l ^. Note that we 
were only able to calculate a nematic-smectic transition 
line for zig-zag molecules. It is plotted in the phase dia¬ 
gram of Fig. [T2] (a) as blue dotted line. 

As already discussed, for the zig-zag molecules On- 
sager theory qualitatively captures the behavior of both 
the isotropic-nematic and the nematic-smectic transi¬ 
tion lines. Both lines are underestimated by the theory, 
but that is hardly surprising. Onsager theory relies on a 
second-order virial expansion of the system’s free energy, 
which is strictly valid only in the dilute regime, while in 
our model orientational and translational order occur at 
relatively high densities. Indeed, it was shown that virial 
coefficients of order higher than two are not negligible in 
hard-needle systems^. In addition, Onsager theory does 
not take into account topological defects that play a sig¬ 
nificant role in the two-dimensional phase behavior and 


phase transitions. In Ref. [13 good agreement between 
predictions from Onsager theory and results from Monte 
Carlo simulations are reported for zig-zag molecules. We 
believe this agreement to be due in part to the small 
system size of 200 molecules, where defects cannot fully 
develop. 


VI. CONCLUSION 


Despite its simplicity, the bent-needle model discussed 
in this paper shows a variety of liquid-crystal phases and 
gives an example of how molecular geometry controls 
their formation. 

Chiral zig-zag molecules assume quasi-nematic and 
smectic phases, depending on density and the central-to- 
tail angle a. We use the orientational correlation function 
52 (r), which decays exponentially in the isotropic phase 
and algebraically in the orientationally ordered phase, 
to identify the isotropic-quasi-nematic transition. In 
two-dimensional systems with short-ranged interactions, 
one expects the transition to take place via disclination 
unbindingiS. Indeed, apart from packing effects, the cor¬ 
relation functions described in Sec. IIV Al behave similarly 
to the ones already reported in other two-dimensional 
anisotropic models with hard-core interaction a^^i^"^~ — . 
For the zig-zag molecule the isotropic-nematic transition 
density piN increases with increasing central-to-tail an¬ 
gle a, while the nematic-smectic transition, pnSi exhibits 
the opposite trend. Such behavior is qualitatively cap¬ 
tured by Onsager theory although both piN and pNS are 
underestimated by the theory, as discussed in Sec. [Vj 

Furthermore, we observe a smectic C phase, where the 
central segment of the zig-zag molecule is tilted against 
the layer normal. The clear formation of smectic layers, 
however, indicates that our systems are too small to ob¬ 
serve the dislocation unbinding scenario predicted in Ref. 


Achiral bow-shaped molecules have a much richer 
phase behavior. It can be divided into three regions, 
depending on the value of a. For small a molecules form 
isotropic, quasi-nematic, and modulated-nematic phases, 
as density increases. The isotropic-quasi-nematic transi¬ 
tion is defect driven, which makes the orientational corre¬ 
lation function switch from an exponential to power-law 
decay with increasing density. No polar order is found 
in the quasi-nematic phase. Further increasing density 
destabilizes the orientational order of the quasi-nematic 
phase. A modulated-nematic phase then takes over, 
wherein bow-shaped molecules form layers of supramolec- 
ular arches. The orientational correlation function 52 )?') 
shows a reentrant exponential decay that corresponds to 
the development of the supramolecular arches. At inter¬ 
mediate a, a direct transition from the isotropic to the 
modulated-nematic phase takes place. The modulated- 
nematic structure becomes less pronounced with increas¬ 
ing a, i.e., the layer thickness decreases towards the 
molecule length, up to the point where the formation 
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of supramolecular arches becomes unfavored and anti- 
ferromorphic domains develop instead. For a > tt/S, the 
curvature radius of the molecules becomes too small to 
induce spontaneous bending of the nematic director and 
quasi-nematic order reenters. Anti-ferromorphic smectic 
order is then found at even higher density. 

A very appealing result of our investigation is the iden¬ 
tification of the modulated-nematic phase made from 
supramolecular layers, which strongly depend on molec¬ 
ular geometry. Our simulations are the first to clearly 
demonstrate such a layered structure in two dimensions 
and to relate it to molecular geometry. We find the opti¬ 
mal value of the curvature radius to be around a ~ tt/IO 
when the modulated-nematic phase occurs at the lowest 
packing density. This finding might be particularly useful 
in developing novel functional optical materials based on 
organic bent-core liquid crystal s where the forma¬ 
tion of polar domains can be used to tune the nonlinear 
optical properties of the material. 

Our investigations further illustrate the richness 
of structures, including supra-molecular organization, 
formed by self-assembling particles of different shapes. 
Being able to control molecular geometry thus offers the 
possibility of designing novel materials, in particular in 
two dimensions, and of tuning their properties accord¬ 
ingly. 
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